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Abstract 



A calculation of the next-to-leading order exclusive extrinsic charm quark 
differential distributions in deeply inelastic electroproduction has recently 
been completed. Using these results we compare the NLO extrinsic contri- 
butions to the charm structure function ^(x, Q 2 , m 2 c ) with the corresponding 
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NLO intrinsic contributions. The results of this analysis are compared with 
the EMC DIS charm quark data and evidence for an intrinsic charm compo- 
nent in the proton is found. 
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1 Introduction 



The European Muon Collaboration (EMC) has played an important role 
in the study of charmed quark production in deep inelastic scattering. Their 
results established that photon-gluon fusion in QCD, the analogue of the 
Bethe-Heitler reaction in QED, explains most of the charmed quark con- 
tribution to the deep inelastic structure function F 2 . The scattering cross 
section for the reaction e~(li) + P(p) — > e~{l 2 ) + c(Pi) + X is conventionally 
written clS cL differential in x and y as 



d 2 a 2na 



2 



dxdy Q 4 



■S 



{l + (1 - y) 2 } F 2 (x, Q\ m\) - y 2 F L (x, Q 2 , m 2 c )} (1.1) 



with q — l\ — l 2 , Q 2 = —q 2 , x = Q 2 /2p ■ q and y — p- q/p ■ 1%. The dependence 
on the charm quark mass, m c , is shown explicitly. The square of the center- 
of-mass energy of the electron-proton system is S. The above formula holds 
after integration over the azimuthal angle between the plane containing the 
incoming and outgoing electron and the plane containing the incoming proton 
and outgoing charmed quark. In a typical experimental extraction of F 2 from 
a, the F L contribution is either neglected or approximated by leading order 
QCD. 

The lowest order (LO) photon-gluon fusion model in QCD is based on 
the twist-two term in the operator product expansion which incorporates 
the factorization theorem || for hard scattering. The EMC F 2 (x,Q 2 ,m 2 ) 
data was not in complete agreement with the predictions of this model at 
all x and Q 2 . The disagreement at large x and Q 2 substantiated the claim 
that a second (higher twist) component of charm production was necessary, 
called intrinsic charm (IC) to distinguish it from the twist-two mechanism, 
referred to as extrinsic charm (EC). 

In the analysis of the EMC data |T[, a relatively simple model was used 
for both the EC and IC components. At the time, only the LO contributions 
to F 2 (x,Q 2 ,m 2 ) from both models were available. Afterwards, Hoffman and 
Moore || calculated the next-to-leading order (NLO) corrections to the IC 
component and discussed their effects on the EMC analysis. Based on LO 
photon-gluon fusion, they found evidence for an 0.3 % IC component in the 
proton. Due to the lack of NLO calculations of the EC component for the 
analysis of the EMC data, the previous results are inconclusive. 
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The NLO corrections to the extrinsic component are now also available || 
H) S B 3 • These NLO calculations yield large cross sections in specific regions 
of phase space since new mechanisms occur, such as the t-channel exchange 
of massless gluons. Hence cross sections for "massive" quark production in 
NLO should really not be compared with the LO prediction via the standard 
K factor, although this is often done in the literature. We will comment on 
this later. 

The single particle inclusive results |7| for the charmed quark distributions 
were applied to the EMC data in |L0| without addressing the IC component. 
However, even the single particle inclusive NLO results are not enough for 
a complete reanalysis of the EMC data since the invariant mass of the cc 
pair, M C c, was included in the scale of the running coupling constant and 
cuts were made on the energies of the decay muons. Thus an exclusive 
NLO calculation, retaining all possible information on both the charmed 
and the anticharmed quarks, is required. The NLO corrections to exclusive 
production were derived in [f|, [|, making it possible to include M c5 in the 
running coupling constant. Since the complete NLO results are now available 
for both EC and IC production, a more detailed QCD analysis of the EMC 
results is finally possible, allowing us to make a more reliable determination 
of the IC content of the proton. 

To make this paper self-contained, we begin with a review of intrinsic 
charm and its NLO contribution to the structure function F 2 (x, Q 2 ,ml) in 
Section 2. We follow this with a brief review of extrinsic charm and its NLO 
contribution to F 2 (x, Q 2 , ml) in Section 3. In Section 4 we discuss the results 
of our analysis and conclude that, even with the addition of the new NLO 
EC contributions and the use of the most recent parton densities, an IC 
component is still needed to fit the EMC data. 



2 The Intrinsic Charm Component 

The QCD wavefunction of a hadron can be represented as a superposition 
of quark and gluon Fock states. For example, at fixed light-cone time, 
t = t + z/c, the proton wavefunction can be expanded as a sum over the com- 
plete basis of free quark and gluon states: \^ p ) = J2 m \ m ) 4>m/ P (xi, kr,i, K) 
where the color-singlet states, \m), represent the fluctuations in the proton 
wavefunction with the Fock components \uud), \uudg), \uudcc), etc. The 
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boost-invariant light-cone wavefunctions, ip m / p (xi, kr,i, A«), needed to com- 
pute probability distributions, are functions of the relative momentum coor- 
dinates Xi = kf / P + and kx,i- Momentum conservation demands Y^l=\ x % — 1 
and Y^i=i kr,i — 0, where n is the number of partons in a Fock state \m). 
When an interaction occurs, the coherence of the Fock components is broken 
and the fluctuations can hadronize, forming new hadronic systems For 



example, intrinsic cc fluctuations || can be liberated provided the system is 
probed during the characteristic time, At = 2p\ a \ ) /M 2 Zl that such fluctuations 
exist. 

There is substantial circumstantial evidence for the existence of intrinsic 



cc states. Leading charm production in irN [12] and hyperon-iV [13] collisions 
requires a charm source beyond leading twist. The leading charm produc- 
tion can be explained by the coalescence of intrinsic charmed quarks in the 
projectile wavefunction with spectator valence quarks |14|, [15]. Final state 



coalescence mechanisms [fL6|] cannot consistently describe data with meson 
and baryon projectiles. The NA3 experiment has shown that the single J/ip 
cross section at large Xp is greater than expected from gg and qq production 



17, nil. This experiment has also measured a significant number of corre- 



lated J/ip pairs |19| which carry a larger fraction of the projectile momentum 
than expected from leading-twist QCD. Intrinsic heavy quark states such as 
\udcccc) can explain the production of fast J/ip pairs [20]. Additionally, in- 
trinsic charm may account for the anomalous longitudinal polarization of the 
J/ip at large xf f2l| seen i n — > J/ipX interactions. 

Microscopically, the intrinsic heavy quark Fock component in the proton 
wavefunction, \uudcc), is generated by virtual interactions such as gg — > QQ 
where the gluons couple to two or more valence quarks. The probability for 
cc fluctuations to exist in a light hadron thus scales as a 2 (m 2 )/m 2 relative 
to leading-twist production [Q. Therefore, this contribution is higher twist, 
suppressed by 0{l/m 2 ) compared to extrinsic production. 

The dominant Fock state configurations are not far off shell and thus 
have minimal invariant mass, M 2 = J2i m Ti/ x i where m T:i = \jk 2 -. i + mf is 



the transverse mass of parton % in the state. Intrinsic cc Fock components 
with minimum invariant mass correspond to configurations with equal ra- 
pidity constituents. Thus, unlike extrinsic heavy quarks generated from a 
single parton, intrinsic heavy quarks carry a larger fraction of the parent 
momentum than the light quarks in the state M. It was shown that large 
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%f virtual cc or lepton pairs can be liberated by a relatively soft interaction 



II]. For soft interactions at momentum scale /i s , the intrinsic heavy quark 



cross section is suppressed by a resolving factor proportional to nl/Am 2 JT3 



for hadroproduction and (fi 2 + Q 2 )/(4:m 2 + Q 2 ) for electroproduction (TT|, ^2 

The general form of the Fock state wavefunction appropriate to any frame 
at fixed light- cone time is 

where T is a vertex function, expected to be a decreasing function of w? h — M 2 1 
a measure of the "off-shellness" of the Fock state fluctuation. The vertex 
function is assumed to be relatively slowly varying; the particle distributions 
are then controlled by the light-cone energy denominator and phase space. 
This form for the higher Fock components is applicable to an arbitrary num- 
ber of light and heavy partons. The Fock states containing charmed quarks 
can be materialized by a soft collision in the target which brings the state 
on shell. In the limit of zero binding energy, \1/ is singular and the fractional 
momenta peak at a?, = m^/m^. Note that the denominator is minimized 
when the heaviest constituents carry the largest fraction of the longitudinal 
momentum. The parton distributions reflect the underlying shape of the 
Fock state wavefunction. Assuming it is sufficient to use a mean value of 
to calculate the x distributions, the probability distribution for a general 
n-particle intrinsic cc Fock state as a function of x is 

(/P " N n [al{M^f tJ^-Fr^i^ , (2-2) 



dxi---dx n s {m 2 h -J2i =1 {mf/xi)) 2 



where m 8 = ym- + (fcf^) is the average transverse mass and N n normalizes 
the Fock state probability. In the heavy quark limit, m c , 77% 3> m^, m q and 
the probability reduces to 

Integration over the light quarks and anticharmed quark momenta for 
n = 5, the minimal intrinsic charm Fock state of the proton, \uudcc), gives 
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the intrinsic charmed quark density distribution as a function of the charmed 
quark momentum fraction, 



c(x) oc 



dPj C {x) 
dx 



-N 5 x 2 [-(1 - x)(l + lOx + x 2 ) + 2x(l + x) In a;] . (2.4) 



If there is a 1% probability for intrinsic charm in the nucleon, as previously 
suggested ||, then N 5 = 36. The charmed quark structure function at leading 
order F% f (x) is given by 



<o), 



x) = 8xc(x)/9 



(2.5) 



in the limit where the charmed quark mass is negligible. In a complete 
analysis of the EMC charm data, the massless result is clearly inapplicable. 

Hoffmann and Moore Q incorporated mass effects into the above analysis. 
They first introduced a mass scaling variable £ = 2ax[l + (1 + Apx 2 ) 1 ^ 2 }^ 1 
where p = m 2 /Q 2 , a = [(1 + 4A) 1 / 2 + l]/2 and A = m 2 c /Q 2 . The proton mass 
is denoted by m p . The cc mass threshold imposes the constraint £ < 7 < 1 
where 7 = 2ax [1 + (1 + Apx 2 ) 1 ^ 2 }^ 1 . Then Q2.5p is replaced by 



(2.6) 



F 2 (0) (x,Q 2 ,m 2 ) = 8£c(£, 7 )/9. 



with c(z, 7) = c(z) — zc(j) /7 for z < 7 and zero otherwise; c(z) is defined in 
( |2.4| ). Hoffmann and Moore found that further mass effects could be incor- 
porated by generalizing the operator-product expansion analysis to include 
both the charmed quark and target masses. The final LO result, c.f. eq. (18) 
in HI, is then 



Fi°\x,Q 2 ,m 2 c ) 



ix 



where 



2px 



1 + 4px 2 ) J 't 



9(l + 4px 2 ) 3 / 2 



;i+4A) 



v 1 



c(£,7) + 3^(£, 7 ) 



1 + 2pxt + 



2\x 



The NLO IC component of the structure function is given by 



9 hh z 



(2.7) 



(2.8) 



(2.9) 
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When the lowest order cross section is normalized to 

a { 2 °\z,X)=5(l-z), (2.10) 
the NLO QCD corrections to the IC contribution are given by eq. (51) in 

4 1 ) ( z , A) = ^5(1 - z) U In A - 2 + VT+4XL + [3L 2 
o7r >• V1 + 4A 

+ AL + AU 2 (-d/a) + 2L\nX -4Lln(l + 4A) + 2U 2 (d 2 / a 2 )]} 



+ 



a. 



A- 



1 



3tt(1 + 4A2 2 ) 2 M1-(1- X)z} 2 
x [(1 - z)(l - 2z - 6z 2 + 8z 4 ) + 6Xz(l - z){3 - 15z - 2z 2 + 8z 3 ) 
+ 4A V(8 - 77z + 65z 2 - 2z 3 ) + 16AV(1 - 21z + I2z 2 ) - 128AV] 
2L 



VI + 4Xz 2 



[(1 + z){\ + 2z 2 ) - 2A^(2 - llz - llz 2 ) - 8AV(1 - 9z)\ 



(1-z). 



where 



k 4 (l + 4A) 2 _ 4^ 4 (1 + 2A)(1 + 4A) 2 L 
Vl + 4A^ 2 (l-z) + 

4Xz[l - (1 - X)z 



}. 



In 



(1 + 2Xz + Vl+lAi 2 ) 2 



(2.11) 



(2.12) 



The leading-logarithmic approximation to this formula, given by eq. (54) in 
the same paper, is 



5? } (*,A) = ^ 

37T 



a. 



(l-2z-6z 2 ) 2(l + z 2 )\nz 



2(1 + z 2 ) In A 



1-z 



a-<) + 2(1 + " 

-5(l-z) [3 lnA + 5 + 27r 2 



ln(l - z) 



(2.13) 



The implementation of the plus distributions in eqs. (|2.11 ) and (|2.13|) is 
standard |2T 



dxf(x) 



1 — x 



1 dx^ ffli + /(l) ln(l - a) (2.14) 

1 — x 
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W)(^) + ^/>^ln (1 -*) 

+ i/(l)ln 2 (l-a). (2.15) 

In our calculations of the IC contributions to F2(x,Q 2 ,m 2 ) we have as- 
sumed a 1% probability for IC in the proton. To compare with the original 
EMC analysis, we take m c — 1.5 GeV and fix the renormalization scale in the 
one loop running coupling constant to /i 2 = Q 2 + 20 GeV 2 and Aq^, d = 0.5 
GeV. Later, when the EC calculations are updated with more recent parton 
densities, we take the same value of J^q C d as m ^ ne parton densities and use 
the two loop running coupling constant to be consistent with the EC results. 

Figure 1 shows the LO and NLO IC contributions to F 2 (x,Q 2 ,m 2 ) for 
two Q 2 values, Q 2 = 7 and 70 GeV 2 . The upper dotted line shows the 
massless (Q 2 independent) result, (2.5), the dot-dashed line shows the £ 
scaling formula, (2.6), and the upper solid line shows the full result, (2.7), 
all at LO. At low Q 2 there is substantial difference between the simple form 
given in (2.5) and the kinematically corrected formula in (2.7). The three 
curves are nearly indistinguishable when Q 2 is large, as seen in Fig. 1(b). 
In the same figure we also show the NLO corrections calculated using eq. 
(2.9). These are the leading-logarithmic approximation eq. (2.13), dotted 
lines, and the complete result, eq. (2.11), dashed lines. These two formulae 
yield different curves at small Q 2 but are almost identical at large Q 2 . 

We note here that we have uncovered an error in the implementation of 



eq. (2.9) in recent work f24" |: the constant term in eq. (2.14), /(l) ln(l — a), 
was omitted. This term is due to the non-zero lower limit of integration 
in eq. (2.9) and has a large numerical effect. Contrary to [pi], our results 
show that the heavy quarks evolve faster than the massless quarks using the 
standard DGLAP equations, as suggested by the original results of Hoffmann 
and Moore [g. 

In the remainder of this paper, for the IC results at NLO, we use the 
sum of the kinematically corrected LO result, eq. (2.7), and the full NLO 
correction obtained by using eq. (2.9) with eq. (2.11). The result is given by 
the lower (at x = 0.4) solid line in Fig. 1. Note that this total is actually 
negative at large x and Q 2 , indicating that even higher order terms are needed 
in the perturbation expansion in this region. 
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3 The Extrinsic Charm Component 



Order a s QCD corrections to the charmed quark structure functions for inclu- 
sive production in deep inelastic scattering were first presented in [|]. Heavy 
quark production was assumed to be extrinsic so that F 2 and Fl could be 
calculated from the inclusive virtual photon-induced reaction j*(q) +P(p) — > 
c {Pi) + X- The parton level interaction is 



y*(q) + a x (A;i) -> c(pi) + c(p 2 ) + a 2 (/c 2 ) 



(3.1) 



where ai and a 2 are massless partons and p 2 and /c 2 are integrated over in 
the inclusive calculation. The structure functions can then be written as 

Q 2 «,(/i 2 ) 



F k (x,Q ,m c 



+ 



A-K 2 m 2 c J$ min £ 



+ E/</p(£,A* 2 



i=q,q 



e l f C fej + C fcJ 111 ^ 2 



+ e ? 4i + e c e* o| 1} 



}. 



(3.2) 



where k = 2,L. The lower boundary on the integration is £ m j n = x(4m 2 + 
Q 2 )/Q 2 - The parton momentum distributions in the proton are denoted 
by fi/p(£,,H 2 ) where //, the mass factorization scale, has been set equal to 
the renormalization scale in the running coupling constant a s . Finally, 

jj. 1 ] and ojj, 1 ] are scale independent parton coefficient 
72T) the coefficient functions are distinguished by their 



and Oft , (/ = 0, 1 
functions. In eq. 

origin: cj^ and originate from the virtual photon-charmed quark coupling 
and therefore appear for both charged and neutral parton-induced reactions; 

arise from the virtual photon-light quark coupling; ojj, 1 ] come from the 
interference between these processes. All charges are in units of e. We have 
included the terms proportional to e c ti in (3.2) even though they do not 
contribute to the total partonic cross section. We have also isolated the 
mass factorization scale dependent terms, proportional to cjj. 1 ] . Finally, note 
that eq. ( |3.2|) only holds for Q 2 > 0. In the photoproduction limit there are 
additional terms involving the parton densities in the photon [J2^, |26], |27], 
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Recently, two of us have reported on the results of a calculation of the 
NLO corrections for heavy quark exclusive distributions at fixed x and Q 2 
[§, [J. This allows us to study correlations between the outgoing particles 
in eq. ( |3.1| ). The exclusive results are needed for our reanalysis of the EMC 
data because in the original analysis |J, the strong coupling constant was 
evaluated at the scale /i 2 = Q 2 + M 2 -. Contrary to what was assumed in 
previous analysis of the IC component, M cE is not a constant but a function 
of a; and Q 2 §. 

Since the NLO calculations were not available at the time of the EMC ex- 
periment, it is clearly interesting to make a reanalysis of their data including 
both the NLO corrections and more recent parton densities. For a complete 
description at NLO, we have included the quark and antiquark contributions 
to charm production as well as the gluon contribution. See || for details. 

Our analysis is done in two steps. We begin with the same parameters 
as used in the EMC analysis JT| to confirm that we can reproduce their EC 
results at LO. Using these same parameters, we then calculate the NLO 
contribution to see what effect it would have had. We have used m c = 1.5 
GeV/c 2 and the one loop running coupling constant at scale fi 2 = Q 2 + M 2 Z 
with Aq^, d = 0.5 GeV. In conjunction with the one loop coupling, the scale 
independent gluon density xg{x) = 3(1 — a;) 5 was also used. We will see that 
the EC NLO calculation yields values of ^(x, Q 2 ,m 2 ) which are above the 
data. This can be corrected by adjusting the normalization of the theoretical 
prediction. 

Next we repeat the analysis with modern parton distributions and the two 
loop running coupling constant evaluated at /i 2 = Q 2 + M 2 5 . The CTEQ3 
MS parton distributions with A { q ] cd = 0.239 GeV |29| are used. We also 
compare with MRS and the most recent GRV |JT| parton distributions 
to investigate the dependence on the parton densities. 



4 Results 

We first show the IC and EC results for F 2 (x, Q 2 , m 2 ) as a function of Q 2 for 
fixed v and study the differences between the NLO and LO results with the 
running coupling scale and gluon distribution used by EMC |]J. We retrieved 
F 2 (x, Q 2 , m 2 c ) in bins of x and Q 2 from the Durham-RAL HEP Database 
[ j3^] . The central x and Q 2 values, x and Q 2 , give the average value of u, 
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V = Q 2 /2m p x. The results corresponding to the three curves in Fig. 13 of 
the EMC paper are presented in Fig. 2 for V = 168, 95, and 53 GeV. Note 
that the v bin widths quoted in Q are 160 < v < 220, 100 < v < 160, and 
60 < v < 100, all in GeV. At the values of Q 2 relevant to the EMC data, 
the mass of the charmed quark is not negligible. Thus the charmed quark 
cannot be treated as massless and absorbed in the parton densities. 

Figure 2 shows the EMC data with our calculations of the LO EC con- 
tribution (dotted lines), the NLO EC contribution (solid lines), the LO IC 
contribution (dot-dashed lines), and NLO IC contribution (dashed lines). 
The LO EC contribution agrees with the EMC result, as expected, since we 
used their gluon density. Clearly the EC component can be reduced to im- 
prove the agreement of the NLO results with the EMC data by changing the 
normalization of the gluon density. This also requires some readjustment in 
the normalization of the charged parton densities to maintain the momentum 
sum rule. In any case, the data cannot be reconciled with the EC prediction 
alone since such a readjustment will not fit the data at large V and Q 2 . 

We next show the EMC data for the structure function F 2 (x,Q 2 ,m 2 ) at 
v = 53 GeV in Fig. 3(a) plotted as a function of x. Also in the same figure 
we add the theory predictions for the LO EC contribution, the NLO EC 
contribution, the LO IC contribution and the NLO IC contribution. Again 
we have used the EMC gluon density so that the LO EC result fits most of 
the data at small x while the IC component is important at larger x. Figures 
3(b) and 3(c) show the corresponding results for v = 95 GeV and 168 GeV 
respectively. It is again true that if we readjust the normalization of the 
gluon density to move the theory curves down we still cannot get a good fit 
to the data using only the EC NLO result. 

Rather than elaborate on refitting the old EMC gluon density, we turn 
to a comparison of modern gluon densities with the data. In Fig. 4 we 
again show F 2 (x, Q 2 , m 2 ) for fixed v versus x with the corresponding theo- 
retical curves for the CTEQ3 parton densities |29| obtained from global fits 
to deep-inelastic, Drell-Yan and direct photon production data. Note that 
the EMC data have not been used in the global analysis to obtain the par- 
ton densities. For consistency m c = 1.5 GeV/c 2 was used in both the EC 
and IC contributions. The two-loop a s with Aq^d = 0.239 GeV and match- 
ing across mass thresholds was used for the EC results which are shown for 
scales /i = /zo/2,/1 = /io and \i = 2/i where /ig = Q 2 + M 2 5 . We also show 
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the NLO EC component calculated at fi — //q using the MRS(G) [[SIJ parton 
densities with Aq^d = 0.255 GeV and the GRV parton densities with 

Aq^ d = 0.200 GeV. All parton densities give essentially the same results as 
there is very little difference between the parton densities in this x and Q 2 
range. To demonstrate this we present average values of x in the v bins for 
all the parton density and scale combinations in Table 1. Finally the IC com- 
ponent is shown using the two loop a s {\i) and the values of Aq CD determined 
by the parton densities. It is clear that the data cannot be fit with only the 
EC contribution, even with the most recent parton densities. 

We now discuss this conclusion in more detail. To fit the data at NLO 
a considerable enhancement in i^x, Q 2 , m 2 ) is needed at large x, significantly 
beyond that shown in Fig. 4. The K factor, Ff LO (x, Q 2 , m 2 c )/F^°(x, Q 2 , m 2 c ), 
would then have to have a strong x dependence. To illustrate the variations 
between the LO and NLO EC results the K factors, calculated with NLO 
parton densities in the numerator and denominator and with the two-loop 
a s , are shown in Fig. 5 as functions of x for the EMC v values. We have also 
computed the K factors with the LO parton densties (available for CTEQ3 
and GRV) and the lowest order coefficient functions with the one-loop a s 
in the denominator and the NLO parton densities with the two-loop a s in 
the numerator, as advocated in [f[0[| . Although some differences in the K 



factors are observable, they are not large enough to affect our conclusions. 
We want to remind the reader that the NLO corrections to heavy quark pro- 
duction contain dynamical mechanisms which lead to large enhancements in 
the cross sections [^, |33], [34]] at NLO. An analogous situation in QED would 
be a comparison of the muon pair production cross section in the reaction 
e + e~ — > fi + fi~ in LO, where there is an s-channel pole, with the NLO muon 
pair production by the two photon mechanism, e + e~ — > e + e~fi + fi~, which 
has a t-channel pole. At large energies, the two photon mechanism clearly 
yields a larger cross section but this does not violate the convergence of the 
perturbative expansion for QED. Therefore it is unwise to place too much 
emphasis on a comparison with the LO results. 

The increase in the K factor appears at large Q 2 and large x, close to 
threshold, as may be expected J7|. However, the most important issue is 
whether or not the K factor is a constant in the region of the EMC data. A 
constant K factor would lead to a uniform enhancement of the EC component 
without improving the fit at large x. Indeed, the K factors change by 10- 
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30% for the higher v bins at large x, not enough to significantly affect the 
results. In the lowest v bin, the K factor increases by nearly a factor of two, 
but the EC results are steeply falling and are negligible compared to the IC 
contribution for x > 0.3. Therefore, in the x region covered by EMC, we see 
that the K factor does not change enough to affect our conclusions regarding 
the need for an IC component. 

In Fig. 6 we show the corresponding K factors for the IC component. 
The opposite behavior to the EC component is observed-the K factor falls 
with x, flattening at large x. This is also clearly a threshold effect since the 
largest K factors appear at the highest V and lowest x. In this region, close to 
threshold, the IC component is negligible compared to the EC contribution, 
thus the large K factor will not affect our conclusions. Where the IC and 
EC components are comparable, the IC K factor is small. Note also that the 
IC K factor becomes less than unity for x > 0.2 depending on i>, as can be 
anticipated from the Fig. 1. 

Now let us try to get a good fit to the data with the sum of an EC and an 
IC component to the structure function. We have performed a least squares 
fit to the data using the Levenberg-Marquardt algorithm |[35|| . The normal- 



ization of both the IC and EC components are taken as free parameters, 

F 2 c (x, Q\ m 2 c ) = a ■ F 2 c ' EC (x, Q 2 , m 2 ) + (3 ■ F 2 c ' IC (x, Q 2 , m 2 ) , (4.1) 

with the scale fi = fiQ. The shift in the normalization of the EC component 
may be considered as an estimate of the size of the NNLO contribution, 
which is equivalent to a shift in the scale /i. Since we have already assumed 
a 1% normalization of the IC component, the fitted (3 is the fraction of this 
normalization. The results are presented in Table 2. The errors quoted in the 
table correspond to a 95% confidence level on the central fit parameters. The 
final results for the combined model of eq. (4.1) are shown in Fig. 7 for the 
CTEQ3, MRS(G) and GRV94 sets of parton densities. The table shows that 
given the quality of the data, no statement can be made about the intrinsic 
charm content for v = 53 and 95GeV. However, for v = 168 GeV an intrinsic 
charm contribution of (0.86 ±0.60)% is indicated. For completeness we have 
examined the influence of a resolving factor, R(Q 2 ) = (/i 2 + Q 2 )/(4m 2 + 
Q 2 ) where /i 2 = 0.2 GeV 2 for the IC contribution. This does not alter our 
conclusions. 

We close with some remarks on other recent results. We have already 



mentioned that another group [10| also checked that recent parton densities 
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fit the EMC data, using slightly older GRV parton densities [|3lJ. They 
concentrated on the x and Q 2 regions where there is clearly no IC component 
in an attempt to use F 2 (x,Q 2 ,m1) for a direct determination of the gluon 
density in the proton |]37|]. Therefore there is no overlap with our work. 
A recent analysis of muon production from IC decays at HERA (and at 
fixed target facilities) suggests that even an 0.1% IC contribution could be 
measurable at HERA p8| . We therefore urge experimentalists to make a 
decisive measurement. 
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Figure Captions 

Fig. 1 (a) The IC contributions to the structure function F 2 (x, Q 2 ,m 2 ) at 
Q 2 = 7 GeV: the massless result (2.5) (upper dotted line), the ^-scaling 
result (2.6) (dot-dashed line) and the kinematically corrected formula 
(2.7) (top solid line). Also shown are the NLO corrections given by 
(2.9), with the leading-log result (2.13) (dashed line), and the full cross 
result (2.11) (lower dotted line). The sum of (2.7) and (2.9) using 
(2.11), represents the total IC contribution to F 2 (x, Q 2 , m 2 ) (lower solid 
line), (b) Same as part (a) for Q 2 = 70 GeV. 

Fig. 2 (a) The EC and IC contributions to the proton structure function 
F 2 (x,Q 2 ,ml) plotted as functions of Q 2 for fixed v = 53 GeV. The 
curves show the predictions from the EC photon-gluon fusion model 
(LO: dotted and NLO: solid lines) and from a 1% IC component ( 
LO: dot-dashed and NLO: dashed lines). The data are from the EMC 
experiment [1[| via Durham- RAL HEP data base []32] . (b) The same as 

(a) for v = 95 GeV. (c) The same as (a) for v = 168 GeV. 

Fig. 3 (a) The EMC data for the structure function F 2 (x, Q 2 , m 2 ) at v = 53 
GeV plotted as a function of x together with the EC and IC LO and 
NLO contributions. The notation for the lines is the same as in Fig. 2. 

(b) Same as (a) for v = 95 GeV. (c) Same as (a) for v = 168 GeV. 

Fig. 4 (a) The EMC data for F 2 (x,Q 2 ,ml) at v = 53 GeV plotted as a 
function of x together with the EC and IC results. The solid lines 
are the NLO EC results with ji = /io/2 (upper), // = /z (middle), 
\x = 2/i (lower). The dotted line is the NLO result for the MRS(G) 
parton densities with /i = /i , while the dashed line is the result for the 
GRV94 parton densities with /i = /i . (b) Same as (a) for v = 95 GeV. 

(c) Same as (a) for v = 168 GeV. 

Fig. 5 (a) The K factors for EC production as a function of x for v = 53 
GeV. The solid lines are for the CTEQ3 parton densities with \i = 2/i 
(upper), /i = /i (middle) and fi = fi/2 (lower). The dotted lines is 
for the MRS(G) parton densities with fi = fi and the dashed line is 
for the GRV94 parton densities with \i = fiQ. (b) The same as (a) for 
v = 95 GeV. (c) The same as (a) for v = 168 GeV. 
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PDF 


n 


z, = 53 GeV 


= 95 GeV 


v = 168 GeV 


CTEQ3 


/io/2 


0.12 


0.12 


0.11 


CTEQ3 


A*o 


0.12 


0.11 


0.10 


CTEQ3 


2/U 


0.12 


0.11 


0.10 


MRS(G) 


A*o 


0.12 


0.12 


0.11 


GRV(94) 




0.13 


0.12 


0.11 



Table 1: 

Fig. 6 The K factors for IC production as a function of x for v = 168 GeV 
(solid line), v = 95 GeV (dashed line) and V = 53 GeV (dotted line). 

Fig. 7 (a) The EMC data for the structure function F 2 (x, Q 2 , ml) at V = 53 
GeV plotted as a function of x together with the fitted results from 
(4.1). The solid line is for the CTEQ3 parton densities, the dotted 
line is for the MRS(G) parton densities and the dashed line is for the 
GRV94 parton densities. The parameters a and (3 are given in Table 
1. (b) Same as (a) for v = 95 GeV. (c) Same as (a) for V = 168 GeV. 

Table Caption 

Table 1 Average x values of the EC component for v bins and the various 
parton densities used in this analysis. 

Table 2 Results of the least squares fit of EC and IC contributions to the 
EMC data according to (4.1). Uncertainties in the fit parameters are 
shown at the 95% confidence level. 
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v = 53 GeV 


v = 95 GeV 


v = 168 GeV 
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a 


P 


a 




a 


P 


CTEQ3 


0.95 ± 0.64 


0.36 ± 0.58 


1.20 ± 0.13 


0.39 ± 0.31 


1.27 ± 0.06 


0.92 ± 0.53 


MRS(G) 


1.02 ± 0.69 


0.34 ± 0.58 


1.38 ± 0.15 


0.32 ± 0.32 


1.47 ± 0.07 


0.79 ± 0.53 


GRV(94) 


1.15 ± 0.77 


0.33 ± 0.58 


1.45 ± 0.16 


0.34 ± 0.31 


1.48 ± 0.08 


0.88 ± 0.53 



Table 2: 
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